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'si" ' Abstract 

^O , 5*3 models offer a low energy approach to describe the observed pattern of masses and 

mixing, of both quarks and leptons. In this work, we first revisit an S3 model with only one 
^^ ' Higgs electroweak doublet, where the flavour symmetry must be broken in order to produce 

^fv . an acceptable pattern of masses and mixing for fermions. Then, we analyse different S3 

models, where the flavour symmetry is preserved as an exact, but hidden symmetry of the 
low energy spectra, after the electroweak symmetry breaking. The latter models require 
the addition of two more Higgs electroweak doublets which are accommodated in an 5*3 
doublet. We also explore the consequences of adding a fourth Higgs electroweak doublet, 
^.' thus occupying all three irreducible representations of 53. We show how the various S^- 

invariant mass matrices of the different models can reproduce the two texture zeroes and 
Nearest Neighbour Interaction matrix forms, which have been found to provide a viable 
and universal treatment of mixing for both quarks and leptons. We also find analytical 
and exact expressions for the CKM matrix of the models in terms of quark mass ratios. 
Finally, we compare the expressions of the CKM matrix of the different S3 models with 
the most up to date values of masses and mixing in the quark sector, via a x^ analysis. 
We find that the analytical expressions we derived reproduce remarkably well the most 
recent experimental data of the CKM matrix, suggesting that 5*3 is a symmetry of the 
quark sector. 



1 Introduction 

The Standard Model (SM) has successfully described the fundamental interactions of elemen- 
tary particles. It has nineteen free parameters, most of which belong to the masses of fermions 
and their mixing. Additionally, the introduction of more parameters becomes necessary when 
the massive nature of neutrinos is considered. 

The observed mass spectrum, the mixing pattern, the fact that there appear to be only three 
generations of matter, the origin of Charge-Parity (CP) violation, among other puzzles, lack 
an explanation within the theory, and are generically known as the flavour and CP problems 
respectively (see for example [HIS])- 

The common approach in attempts to solve the flavour puzzle, was by the addition of a 
horizontal symmetry acting on family space in a non-trivial fashion. The symmetry group which 
relates families in a non-trivial way is known as the family or flavour symmetry group. On the 
other hand, it was also noticed that without adding a family symmetry, one can introduce 
texture zeroes in different positions of the mass matrices to obtain concise relations between 
mixing angles and mass ratios (see refs. [HISHG] for recent reviews). In the late nineteen sixties 
[TIIH] a relation between the Cabibbo angle and a quark mass ratio, 6*^2 ~ \/mrf/ms, was found. 
Then, in the seventies, in a series of papers [914T5]. the importance of this relation was realised 
and generalised in order to relate other mixing angles to mass ratios. We, as other authors, 
combine both approaches as a stepping stone to find analytical expressions for the mixing angles 
as functions of the mass ratios. 

The introduction to the SM of a non-Abelian discrete family symmetry is the simplest way to 
relate families non-trivially. The smallest group among these symmetries is the permutational 
symmetry of three objects, S^. We remind the reader that there are basically two types of 
models based on the group 53. First, there are those models which have only one Higgs field 
which is a doublet under SU{2)l and a singlet under S3 [T61 - I27] . In these models the S3 flavour 
symmetry must be broken in order to produce an acceptable pattern of mixing for fermions. 
Second, there are the class of S'3-invariant models in which the 5*3 symmetry is preserved as an 
exact but hidden symmetry of the low energy spectra after the electroweak symmetry breaking 
[2SH5T|. The latter models require the addition of at least two electroweak doublet Higgs flelds, 
besides the Higgs field of the SM. Of these three SU{2)l doublets, two of them are assigned to 
the doublet irreducible representation (irrep) of S3, and the third one is assigned to the singlet 
one. It is also possible to add extra Abelian discrete symmetries, Z„, to further reduce the 
number of parameters. 

In the last decade and the first years of the present decade, the experimental knowledge 
about the magnitudes of all nine elements of the Cabibbo-Kobayashi-Maskawa (CKM) quark 
mixing matrix, as well as the Jarlskog rephasing invariant, have had a remarkable improvement 
in precision and quality (5^. At the same time, we have witnessed a spectacular improvement 
in the determination of mixing and squared mass differences in the neutrino sector [53|[5^. At 
present, it is crucial for the success of a model of quark and lepton mixing to either agree with 



the experimental information with great accuracy or, better, to predict accurately the observed 
mixing and mass patterns. 

Here, we build various 53 models based on the two aforementioned types, i.e., a model with 
the SM Higgs as an 5*3 singlet and models with a total of three or four Higgs SU{2)l doublets 
assigned to different irreducible representations of 6*3. We then compare the quark sector of 
the models with the most up to date data on quark masses and mixing. Since S3 models do 
not offer an explanation of the value of the quark masses, in order to confront the theoretical 
form of the CKM mixing with the experimental data, we need to perform a x^ fit, where the 
observables should be four independent parameters of the CKM matrix, and the parameters to 
be adjusted should be the quark mass ratios and one free parameter for each type of quarks 
(up and down). The mass ratios are not treated as free parameters, since we allow their values 
to vary within the three sigma region given by the best fit to experimental data measurements 
as given in the PDG [52 | I55]. and our computation at Mz- 

We find a remarkable good quality of the x^ fits for the different models in the allowed 
parameter space of each model. This will make it possible to discriminate among different 
models when measurements in the quark sector further improve. 

In a follow up work |l56j, we will confront the corresponding masses and mixing of the 
leptonic sector of each model with the most up to date experimental data. It is important to 
mention that our approach is at low energies, hence making no assumption of an ultraviolet 
completion of the theory. However, some of these scenarios can be embedded in Grand Unified 
Theories (see for instance refs. |57f[5U] ). 

The present work is organised as follows. In Sections |2]|U we present the basic ingredients 
of the 53 models that we confront with the experimental data on quark masses and mixing. 
In Section [5l we present the form of the quark mass matrices and relate them to two texture 
zeroes or Nearest Neighbour Interaction (NNI) mass matrices. In this way, we are able to 
derive explicit analytical expressions for the elements of the CKM mixing matrix in terms of 
the quark mass ratios and a few free parameters. In Section [6l we present the prediction of 
the CKM matrix for each model. In Section [71 we present a detailed x^ analysis of the CKM 
matrix elements and comment on the very good quality of the fit of our models to the most 
recent experimental data. We conclude in Section [S] with some remarks and an outlook of the 
present work. 

2 53 as a family symmetry group 

5*3 is the symmetry group of permutations of three objects, which can be geometrically repre- 
sented by the different rotations that leave invariant an equilateral triangle. It has six elements, 
the smallest number of elements in non-Abelian discrete groups. It has three irreducible repre- 
sentations (irreps): a doublet 2, and two singlets. Is and 1a, symmetric and anti-symmetric, 
respectively. The Kronecker products of the irreps are: 1^ ® I5 = I5, 1^ (8 1a = Is, 



1a®15 = 1a, l5®2 = 2, 1a®2 = 2, and 2®2 = 1a©15©2. 

The only non-trivial tensor product is that of two doublets, p^ = {pdi,Pd2) and q^ = 
{(1di,(1d2), which contains two singlets, rs and r^, and one doublet, r^ = (r£)i,r/)2), where 

r^ = PDiqol + PD2qD2, rA = PDiqD2 - PD2qDl, 
^D = irDl,rD2) = {PDiqD2 + PD2qDl,PDiqDl -PD2'?D2)- 

So far, the experimental evidence points to the existence of only three generations of quarks 
and leptons [BT], and we will work under this assumption. The SM Lagrangian makes no 
distinction between the different families when no Yukawa interactions are present. We take 
this as a theoretical suggestion of the possible family symmetry relating the three generations 
of matter. In the following section we will discuss first the matter content and then the Higgs 
field content of the different models presented here. 

3 Matter content of S3 

For fermions, when we choose a three dimensional representation of the group 5*3, according 
to the dimension of the fermion mass matrices, we are led to consider two different direct 
irreducible decompositions: 85 = 2 © l^, or, 3^ = 2 © 1^- Then, the assignment of quark 
families to the irreducible 6*3 representations is suggested by the observed mass hierarchy in 
each fermion sector 



rriu : rric : nit ^ 10 -10 : 1, rrid : m^ : rra^ ^10 : 10 : 1 



rUe : rUf, : nir ^ 10 ^ : 10 ^ : 1. (1) 

Hence, we generically assign the first two families, fi{L,R) and fn(L,R), to the doublet represen- 
tation, 2. Then, the third family, fin{L,R), can be chosen to transform as a symmetric singlet 
representation, Ig [5l], or as an anti-symmetric singlet representation, 1a [SlEl]. Here we 
consider the two possibilities and compare the differences that arise among them. The fermions 
in the doublet representation are denoted by 

fliL,R) \ 2, (2) 

JlIiL,R)J 

where /, // or /// represent the family index of a left- or right-handed fermionic field, /(i) or 
/(/J), respectively. Specifically for quarks we have 

fiiiL = {bL,tL) , fiiiR = tn, or fiiiR = bu 

fiL \ ( {uL,dL) \ f fiR \ ( ur \ f fiR \ ( dn 



fllL J V i'^L,SL) J ' \^ fllR J f^u \ CR J ' V fllR J f=d \ ^R ' ' 

in these expressions, {ul, d^) and {cl, Sl) are doublets under SU(2)l, while m/j, cr, dji, and Sr 
are SU(2)l singlets. 



4 Higgs field content 

A state compatible with a SM-like Higgs boson has recently been observed at the LHC [631465] . 
We do not have yet any experimental information about the scalar sector of the SM at higher 
energies, thus it is natural to ponder what are the consequences of having more than one Higgs 
doublet in extensions of the SM (without SUSY). We thus explore here various scenarios with 
different numbers of Higgs SU{2)l doublets. 

There are many possibilities to form S^ invariant scalars with the fermions in the 5*3 repre- 
sentations of Eq. on]) and the Higgs fields assigned as 

Hdw = { tt ~ 2; Hsw ~ Is; Haw ~ 1a- (4) 



^H2W^ 

The cases we consider are the following: 

I. The SM with addition of an 5*3 family symmetry, where a single Higgs SU{2)l doublet is 
a singlet under 5*3 . This model can only explain fermion masses and mixing angles when 
the 5*3 symmetry is broken [T6 |[T7l[22] . 

II. An 5'3-invariant extension of the SM with three Higgs SU{2)l doublets, either as How 
and Hsw or as How and Haw- The choice of the symmetric or anti-symmetric singlet 
depends on the resulting form of the mass matrices we want to generate. Here we consider 
only the invariant scalars that lead to forms of the mass matrices that are able to reproduce 
the measured values of the CKM matrix, namely the Fritzsch two zeroes texture form, 
and the Nearest Neighbour Interaction (NNI) one. 

HI. An S'3-invariant extension of the SM with four Higgs SU{2)l doublets, which are assigned 
to all three irreducible representations of S3: Hdw, Hsw, and Haw- 

In the following section we discuss the cases II and HI. The general discussion of the case I can 
be found in ref. [22], here we present its general features. 

5 Mass matrices 

5.1 Model with one Higgs doublet 

This model has been thoroughly studied in the literature [TB | [22 | [^. for completeness, we briefly 
present it here to compare it with the other models analysed in this work. In this case, the Higgs 
boson of the SM is an SU{2)l doublet and, since it has no flavour, it is in a singlet representation 
of 5*3. When 5*3 is an exact symmety and the field in the S3 singlet representation is assigned to 
the fields in the third generation, then, in a symmetry adapted basis, all elements in the mass 
matrices should vanish except for the element (3,3). Hence, from the structure of the mass 



matrices, only the ferniion that is assigned to the singlet representation of 6*3 acquires mass, 
and the symmetry should be broken in order to give mass to the other families. Realistic Dirac 
ferniion mass matrices could result from the flavour permutational symmetry S^l ® S^f} and 
its spontaneous or explicit breaking according to the chain: 5*31, ® S^r D 5*3"^^ D S2L ® S2R D 
^2^^"^ [IS1[I7]. 

Under an exact S^l ® S3R symmetry, the mass spectrum for either the quark sector (up or 
down quarks) or the leptonic sector (charged leptons or Dirac neutrinos) consists of one massive 
particle in a singlet irreducible representation and a pair of massless particles in a doublet 
irreducible representation of S^l (8> S^ji. Thus, in the electroweak basis, the corresponding 
mass matrices, M^g , are invariant with respect to a permutation of the family (columns) and 
flavour (rows) indices, and all entries in M^g are equal, see Eq. (2.4) in |16) . Once an explicit 
assignment of particles to the irreducible representations of S3 is made, it is convenient to 
make a change of basis from the electroweak basis to a symmetry adapted or hiererchical basis 
by means of the unitary matrix that diagonalises the matrix M^g , 



where 



M;r=utM;ru, (s) 



U = ^ I -V3 1 V2 \ and M*^' = 771^3 I I . (6) 
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In the electroweak basis the masses for the first two families are generated by introducing the 
terms, 

and Mil = -J^ 

y^r ^ .^11 c.^l^ .^.ji c.^l^ ^ / ^ 

(7) 
In the matrix Mj2, Oi and ft are real numbers that parametrize the most general form of a 
matrix invariant under the permutations of the first two columns or rows. In Mji, An and Ai2 
are also real parameters, through which it is possible to construct a complex representation of 
S'3''^^, that allows us to have a CP- violating phase in the mixing matrix [16]. The term Mj2 
breaks the permutational symmetry 5*3^ ® S^r down to S2L ® S2R and mixes the singlet and 
doublet representations of S3, while the term M^^^ transforms as the mixed symmetry term 
of the doublet complex tensorial representation of the 5*3^^^ diagonal subgroup of S^l ® S^r. 
Thus, taking into account all the terms of Eqs. ([S])-(IZD, the mass matrix M^ in a symmetry 
adapted basis takes the form 

M^'=mis\ A* B, Ci \ , i = u,d,l,p^, (8) 






Ai 


A* 


B, C, 





a D, 




where Ai = \Ai\e'^'^\ Bi = — Aj + 5j, and Di = 1 — 6i. From the strong hierarchy of the masses 
of the Dirac fermions, rrii^ » mj2 > ^^n, we expect 1 — 5^ to be very close to unity. 

5.2 Models with three or four Higgs fields 

After the electroweak symmetry breaking, the Higgs SU(2)l doublets acquire real vacuum 
expectation values (vev's), 

wi = {0\Hiw\0), W2 = {0\H2w\0), vs = {0\Hsw\0), and va = {0\Haw\0), (9) 

giving masses to all fermions of the SM. The notation we use for the effective mass Lagrangian 
of these models is 

£q = ~iuL,CL,iL)M''sA Cr - 1(4, s^, 6l)A^^ J Sr \+h.c.. (10) 

In order to make the analysis as general as possible, we first write all possible Yukawa inter- 
actions which arise with the matter and Higgs content of Eqs. ([3]) and (jl]) when assigning the 
third family to different representations. 

Third family in the symmetric singlet representation: For this case, both fniL and 
fiiiR transform as Is 

— Cyj = Y^ {fiiiLflllRHsw) + ~^^2 iflwfi-R + f2wf2R)HsW + 

If— — — — 

9-^3 [{flwH2W + f2wHlw)flR + iflW^lW — f2wH2w)f2R\ + 
~~f^^ (/1H//2R - f2whR)HAW + ~7k^5 iflwHlW + f2wH2w)flIIR + 

-^y/ {findH.wfiR + H2wf2R)\ + /i.e., / = rf, e, (11) 

where Y, , with j = 1---6, are complex Yukawa couplings. When writing the Yukawa La- 
grangian, for up quarks or Dirac neutrinos, the Higgs fields should be replaced by the respective 
conjugate Higgs fields, Hiw — >■ ia2H*yy, i = 1,2. After the electroweak symmetry breaking, 
both up and down quark mass matrices take the generic form 

lV2Yivs + Yiw2 Yiwi + V2Y/vA v^^V 
Mi, = Y^wi - V2Y/va y/2Yivs - Yiw2 V2Y/w2 | , (12) 

V V2Yiwi V2YJw2 2Y/vs 

whose eigenvalues will be denoted as mj, z = 1, 2, 3. 



Third family in the anti-symmetric representation: In this case, both fniL and fniR 
transform as 1a 

- Cyj = ^1 [f niLfiiiRHsw) + ~7k^2 ifiwfiR + f2wf'2R)Hsw + 

2^3 [UlwH2W + f2wHlw)flR + UlV/HlW - f2wH2w)f2R\ + 

If- - 1/- - 

^^4 {flwf2R - f2wf^R)HAW + ^^5 {flwH2W " f2wHlw)flIIR + 

^^6'^ [finLi.Hiwf2R - H2wfiR)] + h.c, f = d, e. (13) 

Third family in mixed representations: fniL transforming as 1a and fnin transforming 
as Is 

- ^Yf = Yi [fniLfiiiRHAw) + ~7k^2 {fiwfiR + f2wf2R)Hsw + 

2^3 [(/lH/-f^2VK + f2wHlw)flR + (/liy-^lH/ " f2wH2w)f2R\ + 
^^4 (/1H//2R ~ f2wflR)HAW + ^^5 (/iVK-^lH^ + f2wH2w)flIIR + 

TT^^^e^ \f11iLiH1wf2R - H2wfiR)] + /i.e., / = d, e , (14) 

and ////L transforming as Is and finR transforming as 1a 

- Cyf = Yl [f IIILfniRHAw) + ~Jk^2 {flwflR + f2wf2R)HsW + 

2^3 [iflW^2W + f2wHlw)flR + {flV/HlW - f2wH2w)f2R] + 

1^_ _ 1/- - 

^^4 iflwf2R — f2wflR)HAW + ^^5 {flW^2W " f2wHlw)flIIR + 

-^^'e^ [7mL(^w/ii? + ^2W^/2i?)] + /^.c, f = d,e. (15) 



We define the notation 

fx{ = V2Yivs, l-ii = Y/w2, fxi = 2Y/vs, 

^Ii = Yiwu fil = V2Y/vA, fii = V2Y/w,, 

f^^^ = V2Yiw2, fii = V2Yjwu fil = V2Yjw2, 
4 = 2Y/va, 



(16) 



which will allow us to express the mass matrices in a concise way, focusing on the number 
of effective parameters entering into each matrix. In the following subsections we present the 
constraints that we require in order to have a successful description of quark masses and their 
mixing for each of the above mentioned cases. 

5.2.1 Case II: Three Higgs fields 

Cases with Hdw and Hsw- When the left- and the right-handed parts of the third fermion 
family are assigned to different singlet representations of 6*3, symmetric or anti-symmetric, and 
are coupled to a Higgs field in the singlet symmetric representation of 6*3, its Yukawa coupling 
vanishes. Hence, this possibility is not feasible for a model of fermion masses precisely because 
the masses of the third fermion family are the largest ones. Due to this fact, if we choose a 
Higgs field in the singlet symmetric representation of 6*3, both left- and right-handed parts of 
the third family of fermions, fin{L,R), must be chosen to transform either as the symmetric 
or as the anti-symmetric singlet representation of 6*3. We can obtain the form of the mass 
matrices by taking the limit of Haw — )■ in Eqs. (fTT]) -f fT3|) . for the cases of either both left- 
and right-handed parts of //// in the symmetric, or anti-symmetric singlet representation of S3, 
respectively. 

The first two cases, A and A', of Tab. [H corresponds to the case where both the left- and 
right-handed parts of the third family are in the symmetric singlet representation of 5*3, Is- 
The cases B and B', correspond to the case where both the left- and right-handed parts of the 
third family are in the anti-symmetric singlet representation of 6*3, 1^. The first column of 
Tab. [T] denotes the labelling we use, while the third one, gives the form of the mass matrices 
after the electroweak symmetry breaking (EWSB). Note that for these matrices, the elements 
(1, 1), (1, 3), and (3, 1) are not different from zero. The fourth column, sub-case A, corresponds 
to a matrix Aijji^j., where we have rotated the matrix to a basis where the entries (1, 3) and 
(3,1) vanish and we have subtracted the element A^/^jerl-'-' -'-]' which will be denoted by /Iq, 
from the diagonal 

/iQ a^ 

Mi^Mi,,^^n{e\,Min{e)l={ af* bf cf | = /i^isxs + a^L. , (i7) 

cf* df 

where a^ , b^ , c^ , d^ , and Hq are shorthand for the entries in the rotated mass matrix Aig . 
Now, the matrix A^^^jer ^^^ ^^^ texture zeroes 

if 

-A^Le. = I «^* b'f cf 
cf* d'f 






af \ 


af* 


bf -fil cf 





cf* df-4J 



and eigenvalues denoted as a/, i = 1,2,3. Then, the physical masses mi are related to the 
shifted masses af simply by 

^i = /"o + ^I- (19) 

Both transformations, the shift and the rotation, are unobservable in the quark mixing sector, 
as long as we rotate both matrices, in the u and d sectors, with the same angle 9 B The rotation 
in the right hand side of Eq. flTOj) is absorbed in the redefinition of the right-handed quarks 
of Eq. (do]). 

The vanishing of the entries (1,3) and (3,1) in the rotated mass matrix Aif^i^^ i^ only 
possible if the rotation angle, 6, and the real expectation values of the Higgs fields in the 
doublet irrep, Wi and W2, see Eq.Q, are related by the condition 

t&ne = wi/w2. (20) 

The diagonalising matrices that enter in the definition of the quark mixing matrix, Vckm, 
may be obtained from the diagonalisation of ^4^1^^ instead of M.g.^. Since there are no right- 
handed currents in these models, Aif^i^^ '^^Y ^^ constrained to be Hermitian without any loss 
of generality [1]. 

It is interesting to notice that, in order to reproduce the Nearest Neighbour Interaction 
(NNI) mass matrix form [HZIES], it is enough to fix the rotation angle in Eq. (I17p at ^ = 7r/6, 
no Hermiticity of the mass matrix is required. This is also important since the NNI mass matrix 
form has been shown to provide a good description of the mixing angles in a unified treatment 
for quarks and leptons |22|I69]. This is the subcase A' in Table 1. 

When the third family of fermions is assigned to the anti-symmetric singlet representation 
of 5*3 we obtain fermion mass matrices with two texture zeroes or the equally successful NNI 
form. These are the cases B and B' in Table 1. 

This shows that the requirement of invariance under the 6*3 flavour symmetry group gen- 
erates the phenomenologically successful Fritzsch-like mass matrices with two texture zeroes 
form and the equally successful NNI form. 

Cases with Hdw and H^w. As mentioned before, when the left- and right-handed parts of 
the third generation fermion fields are assigned to the singlet representations of S3, one to the 
symmetric and the other to the anti-symmetric one, the Higgs field in the singlet representation 
of 5*3 should be assigned to the anti-symmetric singlet of S3 to form a non-vanishing Yukawa 
coupling. 

The form of the resulting mass matrices is shown in the third column of Tab. HJ The form 
that these matrices take after a transformation to a basis where some of its elements are zero 
is shown in the fourth column of the same Table. The particular cases when the rotation angle 
is ir/Q, denoted as C and D' correspond to the second and fourth row of Tab. [2j 



^The details of the shift and the rotation are given in Section [5T3l and in the Appendix lAl 
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3HDM: GsM ® ^3 



J-'l J-r Mass matrix (FB) 



Possible mass textures 



A 2,1s 



2, 1? 
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,/^2 



(3-t^) 
2/i^c2 (1 -St^) 
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■m{ 







■M2 



c2(l-3i2). 
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v^'^2 U J3M7 

75/^9 Ms - Ml/ 



s 



2, 1/ 



2, 1/ 




m{ M7 

If f 

Ms M3 y 



-m{c2 (1 . 



-Mic2 (1 - 3*2) 



■3*2) 2^|sc(3-t2) 
-Me 7c 



M3^- 





-M^/c 
" lMisc(3- 



*')y 



B 




Table 1: Mass matrices in 6*3 flavour models with three Higgs SU{2)l doublets: Hiw and H2W; which occupy 
the 5*3 reducible representation 2, &ndHsw, which transforms as Is for the cases when both the left- and right- 
handed fermion fields are in the same assignment. The mass matrices shown here follow a normal ordering of 
their mass eigenvalues (to{ , rrij , ^3 ) . We have denoted s = sin0, c = cos0 and t = tan0. The third column of 
this table corresponds to the general case, while the fourth column to a case where we have rotated the matrix 
to a basis where the elements (1, 1), (1, 3) and (3, 1) vanish. The primed cases, A' or B', are particular cases of 
the unprimed ones, A or B, with 9 = 7r/6 ov 9 = 7r/3, respectively. 
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3HDM: GsM ® ^3 



Name J-^ J^r Mass matrix (FB) Possible mass textures 

Ai2 



C 2, 1a 2, U 




112 

^21 ^22 ^23 

-A32 AssCm^-^i'I) 

B12 -B23 

B21 S22 

B32 B33{l^{^ul) 



a 




D 2, Is 2, 1, 



■ 4 


4 + M^' 

-M2 

4 


-m/ 

^3 


, M^ 



D' 





^32 ^33(^3 ->i^:{)y 

A12 A23 \ 

A22 

A32 A33{pl^u[)J 

B12 ^ 

B22 B23 

Bssipl^uf,)^ 



A'.,^ A'^.,ii,f-^uf)j 



K-Al^i ^ "3) 



Table 2: Mass matrices in S'3-invariant family models with three Higgs SU{2)l doublets: Hmr and H2W are 
assigned to the S3 doublet irreducible representation 2, and Haw, transforms as 1a- We consider the cases 
of having the left and right parts of the fermion fields assigned one to the symmetric and the other to the 
anti-symmetric singlet representations of 5*3, since only these combinations give rise to a non- vanishing (3,3) 
entry in the mass matrix. In all cases it is assumed that the mass eigenvalues follow a normal order in terms 
of magnitudes. To simplify the notation, Xy denotes the entry {i,j) in X which is any of the mass textures 
shown in Table [T] The primed matrix elements, ^J- or iJ^' . are particular cases of the unprimed ones, Aij or 



Bij^ with = 7r/6 or = 7r/3, respectively. 
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5.2.2 Effective number of parameters 

Pakvasa and Sugawara j2^ analysed for tlie first time tlie Higgs potential involving two Higgs 
fields in the doublet irrep of 5*3 and a third one in the symmetric singlet irrep. They found 
an accidental 5*2 symmetry at the minimum if one requires wi = W2, which in turn implies the 
following equalities 

/^2 = /^4. /^6 = /^7' d = f^9 ^ (21) 

reducing the number of parameters. The corresponding mass matrices thus reduce to the cases 
of Tab. [3l where we have written down the effective number of free parameters involved in 
each sector. The form of the matrices of the fourth columns of Tabs. [T]l2] is independent of this 
assumption. Here, by effective free parameters we mean the independent real parameters to be 
adjusted on a x^ analysis. Hence, the counting is performed by considering that each matrix 
has only a single independent phase plus the number of absolute magnitudes of the complex 
free parameters. Comparing the number of effective free parameters with the number of real 
positive parameters in the general case. 



/^i: /^2: /^3' ll^rli 6^ and arg /i^ , (22) 

we find they are equal, since in both cases since in both cases the submatrices formed by the 
elements (1,1), (1,2), (2,1), and (2,2) are parameterized by only two different parameters. 
Since these parameters are physically irrelevant, we conclude that from the point of view of the 
quark mass matrices, the assumption Wi = W2 yields the same result as assuming that Wi and 
W2 are related through Eq. fl2U|) . In the cases A(A') and B(B') of Tab. [2] we notice that after 

Name Mass matrix form No. of effective real free parameters 

A 



B 




Table 3: Form of the mass matrices for the cases of Tab. [T] where we have assumed wi = W2- This corresponds 
to the case where the vacuum of the spontaneous symmetry breaking of the EWxS'3 theory has an accidental 
5*2 symmetry. 

reparametrizing the mass matrices in terms of their eigenvalues, the resulting expressions for 
the entries in the CKM mixing matrix expressed in terms of the quark masses are the same. 
Therefore, the quark mixing is insensitive to whether the assignment of the third family is done 
to the symmetric or anti-symmetric singlet. Both cases will lead to the same result. 
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5.2.3 Case III: Four Higgs fields 




In case III all terms of Eqs. flTT|) . f lT3lfT5|) should be present. As in the previous subsection, we 
have identified the transformations that yield (a) a Hermitian matrix with vanishing elements 
(1,3), (3,1) and (1,1), and (b) a NNI matrix form. Now, taking into account all terms in 
Eqs. (ITTt [T8|) . all the possible forms of the mass matrices that we obtain appear in Tab. HJ 

5.3 Diagonalisation Procedure 

We proceed as in ref. [22], where a general matrix with two texture zeroes, representing mass 
matrices of the basic S^ models, 



(23) 



was diagonalised. However, the models of cases II and III have a non-zero entry in the position 
(1, 1). In order to take these models to the form of Eq. (23), we just make a simple shift as 
follows: 

-^Lr. = /^0l3x3 + -ML.- (24) 

J" 

As a starting step in the diagonalisation of the matrices A^)|/jer' ^^ write the above shown 
Hermitian matrix in polar form in terms of a real symmetric matrix Ai^j^^^ and a diagonal 
matrix of phases 

Pf = diag[l, e*"^!^ e^('^^^+'^2^)], (25) 

Mf / A-f 0\ / \Af\ \ 

-ML. - VJ^^^^^Vf = V} Af* Bf Cf]Vf={ \Af\ Bf \Cf\ , (26) 

^3 \ Cf* Df) \ \Cf\ Df ) 

where the phase cpif is fixed by 0ij = arctan(|/i5|/|/i(|), and the phase 02/? remains a real 
free parameter. Then, as usual, the mass matrix A^j^jgr ^^^y ^^ brought to a diagonal form by 
means of an orthogonal transformation, 

-ML. = 0/diag[?(,-?2^1]0j, (27) 

where 5/ = cr^ /a^ are the corresponding real eigenmasses of Aifji^j, and O/ is a real orthogonal 
matrix. Hence, our unitary matrix, which takes us from the hierarchical basis to the basis 
where the matrix M-jji^j, is diagonal, is 



U/ = OjVf. {2i 



14 



Name Tl J^i 



R 



4HDM: GsM ® ^3 



Mass matrix (FB) 



Possible mass textures 



A 2, Is 2, U 



' Ml + M2 




f f 

Mi +M5 

/ / 
Ml -M2 

Mg 



/ M^'^c(3-t2) 

+m{ 



\ 



-Ms 




M7*/c 



M7/C 



/ 



^M2 Ms 



If f f 

75/^8 ^^■i-^^{) 



B 2, 1a 2, h 



B' 






^/.{c2 (1 - 3*2) 





-A«{c2 (1 - 3*2) 

+d 
2 flics (3 - t2) 



-Ms* A 



-2m{ 



M3 



-Me/c 



~2f,l +pi{ 



^■-/.{+/4sc(3-t2)y 



-2m^ 



-2^6, 
Ms - M{y 



C 2, 1a 2, 1. 



4 Ms 
Ms 


mI-hm^ 

m{ - M^ 
-d 



I A12 

^21 ^22 
VA32 





A23 

^33 (Ms -* 



4)> 





B21 





B12 
B22 

Bs2 



B2S 



S33(M3 ^'^s) 



c" 



^12 





^33 (m;^ 



4)) 



"21 




B, 



R,, 



s: 



23 




^3(^3 



'-:(); 



D 2, Is 2, ly 




M4 -t-Ms 

m{ -M2 

Mg 



f A12 A23 ^ 

A21 yl22 

^ yls2 AssCms -^^i)) 



B12 ^ 

S22 ^23 

BssCm^-^i'^'V 



D' 



^21 




^12 ^23 \ I " ^12 

, B21 

42 4s(M3^;^3)/ V^32 



^33 (M3 



^^3)^ 



Table 4: Mass matrices of S'3-invariant flavour models with four Higgs SU{2)l doublets which are assigned to 
all the 5*3 irreducible representations, HAW^y-sw^'^DW- We denote /i{ = v^Fj ws 



4 = Yiw2,^li = 2Y/vs, 



■^ = Yiwi, 4 = y^YivA, ^^ = ^Yiwi, n{ = V2Y/w2, 4 = V^Yiwi, 4 = y^Y^fwi, and ^ - ^y/va- 



In all cases it is assumed that the mass eigenvalues follow a normal order in terms of magnitudes. To simplify 
the notation, Xy denotes the entry {i,j) in X which is any of the mass textures shown in Table 1. We have 
denoted c = cos 6* and s = sin 6'. The only choices of 53 assignments that may produce a viable model are those 
where both left- and right-handed parts share the same assignment, as the first and second cases. The primed 
matrix elements, A^, or B^' , are particular cases of the unprimed ones, Aij or Bij, with 9 ~ tt/Q or 6 ^ 7r/3, 
respectively. 
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We follow the procedure of ref. [22]. Using the three invariants of the generic real mass matrix 

^Hier 



m{ 



Tr[Mk, 



/ 






~f~f 

-0-10-2 . 



1 



mr + inr + 1 



'HierJ i — V"i; "rV"2; "T J- , (29) 

its parameters, \A^, B^ , \C-^\, and D^ may be expressed in terms of the eigenvalues, a/. 



\Af\ 



~f~f 



|2 



ic^r 



l-{Df- i^i) 



D^-K 



rJ 



i^Df-lJ,-a{){Df -fii-ai) 



(30) 



where we have defined /Iq = fiycrl. To simplify the notation, we define the free parameter Sf 
through the following relation 



6 



f 



1-{D^ 



/^oj 



(31) 



which indeed, together with Eq. fl30l) . allows us to write the mass matrix Aig in terms of its 
invariants and just one free parameter 6f 

/ 



^mer 







V 



l-Sf 





a. 



~f~f 

1-5, 



a( + 5f 







^f cfcf 

f 



1-5,^1^2 



in this expression, we have made the following identificationqj 



(32) 



^1 



f 



a 



f 



6 



/' 



^l 



al-6 



f 



(33) 



such that, 6f is a measure of the splitting of the two small masses in the first two families in 
the 5*3 doublet as a fraction of the mass of the third family in the 6*3 singlet. Therefore, the 
following hierarchy among the 6's for the different kinds of fermions: 



1» 6^> 5i> 5d> 5„ 



(34) 



^In order to make a direct comparison with the notation used in previous pubhcations |161I221[^ . a change 
of labels / <->■ i and ^ -H- fi must be done, everything else remains the same. 
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is to be expected. Note that the form of the matrix in Eq. (l32l) is completely analogous to the 
mass matrix discussed in Eq. (17) of ref. [22], just with the replacement 0"j — )■ rrii. Therefore 
the diagonalising procedure will follow exactly as in ref. |22], and consequently, the form of 
the CKM matrix will be the same, we just need to replace rrii by ai and take into account 
the appearance of a new phase (pi. We should bear in mind that a^ are shifted masses, so 
in Eq. (132|) there are three physical invariants involved and two free parameters, Sf and JIq. 
The CKM matrix should contain only one physical phase, the CP violating phase, which means 
that if there are two parametric phases, 0i and 02, the CP violating phase will be a combination 
of both. 

In what follows we describe the general procedure to find the diagonalising matrices for all 
the cases considered. We then proceed to give the specific details for the cases I through III, 
mentioned in section HI 

5.3.1 Case I: A single Higgs field 

In this case, we have only one Higgs field transforming as Ig. There is no Higgs field assigned 
to the anti-symmetric representation, 1a, which translates into the vanishing of (pif. Also, in 
this case there are no shift parameters /ig for f = u, d. 

5.3.2 Case II: Three Higgs fields 

In this case, when the left- and right-handed parts of the fermionic fields of the third family are 
assigned to the singlet symmetric representation of S^, there cannot be a Higgs field transform- 
ing as 1a, therefore the phase (pif vanishes. However, there are shifts /Iq, which in principle are 
non- vanishing. 

5.3.3 Case III: Four Higgs fields 

In this case we have a fourth Higgs field assigned to the anti-symmetric singlet irrep 1a, which 
produces the /ig parameter in the mass matrices shown in Table 4. In cases A and B in Table 
4, we find that after reparametrizing the corresponding mass matrices in terms of the mass 
eigenvalues, the resulting reparametrized mass matrices are equal. In the reparametrized form, 
the following inequality holds A^{f-gj,[2, 2] > A4;^jg^[l,2] and this inequality implies that either 
/ig vanishes or the relation wf = 3u7| is satisfied. 

In this work we will avoid taking a particular value for the rotation angle, and in consequence 
we will assume that /ig vanishes. Therefore, the Higgs transforming as the anti-symmetric 
singlet representation, 1a, does not contribute to the Yukawa matrix and the phase (pif does 
not appear. Hence, the mass matrix, A^^^jg^, has only one CP violating phase 02/ and the 
parameter 6f is now constrained to satisfy 

Gf{6f, a/)t\3 - ef + Aa{ai{l - 31^ = 0, (35) 
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or 



Gf{6f, a/){l - 3t2)2 + 4a(a(t\3 - ^f = 0, (36) 

for cases A or B, respectively, where Gf{Sf, a/) = 5^ — [1 — 2(?( ~5^2)]'^/ + (^i ~^2)(^i ^^i ~ 
2)5 f - ia( - aif. 

6 Form of the CKM matrix 

The VcKM matrix is defined as 

VS'km = Ul,U,, = 0^p(--'^)0,, (37) 

where p("~'=') = diag[l, e'^^^ e«{0i+02)j -^jth 0j = 0^^ — 0^^, and 0^,^ are the real orthogonal ma- 
trices, Eq. (l28l) . that diagonalise the real symmetric mass matrix of Eq. (132|) . The substitution 
of the expressions O/ g in the unitary matrices of Eq. (1371) allows us to express the entries in 



^This is completely analogous to the expression of Eq. (25) in ref. [22], with the replacements rrii — ?> ai and 






the quark mixing matrix V^^a/ ^^ explicit functions of the quark masses 

< = -\l'-iM + \/5fc (v/(i-'^«)(i-'5.)erei + ^^;Aefe?e^^^) e^ 



V 1'3u^3d I V T^-iuT^'id V 'D-iuD'id j ' 



^.. 



with 



^2{«,d) = (1 - 5«,d)(?«,d + ?c,s)(l +5c,s), 

^3(n,d) = (1 - 5«,d)(l - 5'«,d)(l + ?c,.)- (39) 

6.1 Case I: A single Higgs field 

For this case, the form of the CKM matrix corresponds to that of Eq. (|38|1 with 02 = 0, /Xq 
vanishing and m[ = a/. 
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6.2 Case II: Three Higgs fields 

For the cases A and B in Tab. [H we can parameterize the CKM matrix with a non-vanishing 
phase 02 = 0, and since for these cases /Xg is not zero, we set mf = /Xq + aj . 

6.3 Case III: Four Higgs fields 

In this case, we can parameterize the CKM matrix as in cases A and B, again with a non- 
vanishing phase 02 = 0, following the discussion in Sec. (I5.3.3p . In these cases fxl does not 

f f f 

vanish and we set ?7i^ = ^q + af . 

7 The 5*3 models compared with experimental data 

7.1 Experimental status 

Over the last decade, there has been a remarkable improvement in the precision and quality of 
the measurements of the elements of the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix, 
the quark masses, and their uncertainties. At present, any model for quark masses must provide 
a detailed analysis of their predictions or have the ability to reproduce or to accommodate the 
ever increasing precision of the experimental results. In the following subsections, we present 
a brief overview of the experimental status of quark masses and their mixing. We explain how 
we confront this information with the exact analytical expressions we found, given in terms of 
quark mass ratios, of the CKM mixing matrix elements. Then, we comment on the results of 
the x^ fits for the 5*3 models presented in this work. 

In order to do this analysis, we will use the running quark masses at the electroweak scale 
which we fix at the Mz scale. We obtained the numerical values of the running quark masses at 
Mz using the RunDec program (70] and the most recent PDG values [52], which are presented 
at different scales. For comparison, we also quote the results on values of the quark masses 
2010-2011, reported in ref. |S5] and in the 2011 online version of the PDG. We present the 
results in Tab. [5l We note that while most part of the central values of mi{Mz) are compatible 
with those cited in ref. |71| . the uncertainties have been greatly reduced in the last analysis. 

7.2 Fitting procedure 

We construct the x^ function as 

2_ (ld-|K.|)' , (iei-|K.|)' , (ICI-IK.!)' , {Jt-J,T ,,.. 

Vud Vua Vub Jq 

where the quantities with super-index "th" are the complete expressions for the CKM elements, 
as given by the 5*3 models, and those without, the experimental quantities along with their 
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2011 values [GeV] 


2012 values [GeV] 


mf^{Mz) 2011 
[GeV] 


my^{Mz) 2012 
[GeV] 


mt 


172.0 ±0.6 ±0.9 


172.85 ±0.71 ±0.85 


171.13 ±1.19 


171.07 ±1.21 


mos 


4 67+°-^'^ 


4.65 ±0.03 






rrih 


4 10+0.18 


4.18 ±0.03 


2.84 ± 0.04 


2.85 ±0.04 


rric 


-■-•^^-0.11 


1.275 ±0.0025 


0.616 ±0.064 


0.626 ± 0.0094 


nis 


n 1 nn+0030 

U.1UU„Q020 


0.095 ±0.005 


0.061 ±0.015 


0.055 ± 0.0033 


rrid 


(0.0041,0.0057) 


4:.8toi X 10"3 


0.00284 ±0.00050 


0.0028 ± 0.0005 


rriu 


(0.0017,0.0031) 


2.3+J^:^ X 10-3 


0.00139 ±0.00042 


0.0014 ±0.0005 



Table 5: Values of quark masses, in GeV, as appear in the online 2011 version of the PDG and in [55], and 
the updated values of 2012 [52]. The values quoted at M^, were obtained with the program RunDec [70] at 
four loops in the running of as- For the 2011 data, note that with the use of the chiral perturbation relation 
involving the parameter Q = 23 ± 2. we obtain mu(Mz) ~ 0.00130 ± 0.00047, which is now compatible with the 
value obtained directly from the PDG and the evolution up to Mz- 

uncertainty ay- . We consider the following experimental CKM values 



2011: iKdl 

I 'usl 

\Vub\ — ^.^^^^. -.- ^.^ww^t:, 

J = (2.91 ±0.155) X 10-^ 



0.97428 ±0.00015, 
0.2253 ±0.007, 
0.00347 ±0.00014, 



2012 : iKdl = 0.97427 ± 0.00015, 
IK.I =0.2253 ±0.007, 
\Vub\ = 0.00351 ±0.00015, 
J= (2.96 ±0.18) X 10-5 , 



(41) 



which correspond to a unitary CKM matrix in the case of three generations of quarks. Since 
unitary of the CKM mixing matrix is assumed, there is no need to make the fit to the entire 
matrix but only to four observables. The theoretical expressions of the CKM elements are given 
in terms of the mass ratios, m^, Eq. (133|) . or the parameters Jj, Eq. (fT9|) . hence the minimisation 
of the defined x^ is a function of the parameters nii, Su, Sd and cos0i. This means that, as 
a result of the minimisation, there is a best fit value for each of those quantities, for which 
X^ takes the minimum value. The mass ratios rhi are not free parameters. The limits we set 
correspond to their allowed 3a regions. We also test if there is convergence when using the 
(2cr) regions. In the fitting procedure, we used MINUIT from ROOT [72] for the numerical 
minimisation. The values used for the parameters m^ are given in Table 





2011 


2012 


m„ (Mz) 


0.0000082 ± 0.0000027 


0.0000083 ± 0.0000030 


mc {Mz) 


0.0036 ±0.0004 


0.0037 ±0.00008 


rud {Mz) 


0.00098 ±0.00018 


0.00098 ±0.00017 


rris {Mz) 


0.0205 ±0.0056 


0.0190 ±0.0014 



Table 6: Comparison of the values of the mass ratios, at Mz, in 2011 and 2012. 
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mf'{Mz) 






rris 
rriu 


0.059 ±0.0066 
0.0028 ±0.0005 
0.0013 ±0.0005 


m, (Mz) 
rrid [Mz] 
mu {Mz) 


0.0205 ± 0.0026 

0.00098 ±0.00017 

0.0000078 ± 0.0000030 



Table 7: Changes in the masses of the Hghtest quarks when using the average of the theoretical determinations 
of TOs, that is, not including attice determinations. We obtain ms{2GcW) = 0.101 ± 0.011 GeV. 

7.3 Results 

We have proceeded with the niinimisation of the x^ as follows. We used MINUIT and varied 
all the parameters rhi, within the 2cr and 3a ranges given in Tab. [6] and Tab. [71 and 6u, Sd as 
true free parameters. 



7.3.1 Case I: A single Higgs field 

This case corresponds to the well known case of broken S^l ® S^r, in the presence of one singlet 
5*3 Higgs field, which gives rise to an effective mass matrix of the form of Eq. ([S]). Hence, the 
CKM matrix we fit is that of Eq. fl38|) with 02 = and ctj = rhi. 

In one set of fits, we fixed (pi to 7r/2 following previous fits to the quark mass ratios where 
this value was shown to be the preferred one [T6t[T7]. However, when allowing the phase 0i to 
vary in the region cos0i G (0, 1), the quality of the fits is better for a larger value of cos0i, 
than for a small value of 0i. Therefore, we present three sets of fits, one when 0i is fixed to 
7r/2, another when we allow cos0i to vary in the region (—0.5, 0.5) and a third one, where we 
allow cos 01 to vary in the region (0.5, 1.0). 

We recall that minimisations with MlNUlT rely on setting a starting value for the param- 
eters to fit with a seed close enough to the minimum, therefore if the range of variation of a 
particular parameter is large, it is difficult to find its best fit point. Additionally, to check the 
consistency of a minimum, one should remove the limits of the parameters to fit. Unfortu- 
nately, if we perform the fit leaving the parameters m^ completely free to vary without limits, 
the quality of the fits does not really improve, and most importantly, for these cases it turns 
out that the best fit point of m^ is of order 10~^. We mention that the reported values in 2011, 
by the PDG, of rriu and nid quoted an uncertainty of about 30% of the central value, so they 
were difficult to fit. The situation in 2012, particularly for rris, changed since the uncertainty 
in the lattice determinations of rris was reduced down to 5%. In contrast, the theoretical de- 
terminations of rris have an uncertainty of almost 10%. Consequently, in order to assess the 
impact of lattice and theoretical determinations, we also make fits using only an average of the 
theoretical determinations. 
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Case when (p^ = -nil fixed. We find that for the reported 2012 experimental values of 
quark masses and mixing, when allowing the mass ratios to vary within their 3o" ranges, for 
m„, mc, and m^, the x^ function attains a minimum within the corresponding Sex region of each 
of the parameters above. On the other hand, the BFP of m^ lies within its la region. The 
results of this fit are shown in Figs. [T]-[2]and the values of the BFPs are given in Tab. |H1 

For comparison, we have also performed the fit with the experimental results reported in 
2011 by the PDG and with the 2012 results taking into account the theoretical determinations 
of m^. Tab. [71 As we can see from the plots in Figs. [T]|2] for the former set of data, the x^ 
functions does not really attain a minimum as a function of m„ and rric, when they vary within 
their 3o" region. On the other hand, the best fit point of m^ lies within its three sigma region, 
while the best fit point of mg within its two sigma region. When we take into account the data 
from 2012 only with the average of the theoretical determination of m^, we notice that the best 
fit points of vdu lie within their corresponding 3cr region, while the best fit points of m^ and 
rrij, lie within their corresponding 2cr region. The results of this fit are shown in FigsJSHU the 
values of the best fit point are given in Tab. [Si 

Varying cos 0^ in (—0.5,0.5). When allowing the mass ratios to vary within their 3cr 
ranges, we find that the x^ function does not really reach a minimum as a function of m^ and 
rhc^ while as a function of m^ and ms, it does reach a minimum within their corresponding \o 
region. The results of this fit are shown in Figs. [3M1 the values of the best fit point are given 
in Tab. [8] For the fit, when taking into account the average of the theoretical determination 
of rris, we find that x^ attain a minimum only as a function of m^, whose best fit point lies 
within its \o region. The fit, when using the 2011 data, shows that the x^ function attains a 
minimum as a function of m„, m^ and rhf.. For m„, its BFP lies within its Sa region while for 
mrf and mg, their best fit points lie within their corresponding 2a region. 

Varying cos (p^ in (0.5, 1.0). For this fit, when allowing the mass ratios to vary within their 
3cr ranges, we find that the x^ function attains a minimum as a function of all the mass ratios. 
For m^, rhc and m^^, the minimum lies within their 3o" range. For m^ the minimum lies within 
its \a range. The results of these fits are shown in Figs. ^^ the values of the best fit points are 
given in Tab. [HI When considering the average of the theoretical determination of m^, we find 
that x^ attains a minimum only as a function of rhg^ whose BFP lies within its la region. For 
the data of 2011, x^ attains a minimus for rhu and m^ within their corresponding 3a region, 
while rhd within its la region and m^ within its corresponding 2a region. 

General comments on the quality of the fits. From Figs. [T][4|we can see that the ratio 
rhc is not greatly affected by the change in the value of 0i. However, the minimum of x^ as a 
function function of rhc when cos(/)i ~ 0.5, seems to be better behaved as that of a fixed value 
of 01 equal to 7r/2. When we allow rha to vary in the range determined by the uncertainty in 
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Figure 1: Results of the x^ fit as a function of to„ and ttIs for the case where the phase (j>i is fixed to tt/2. 
The reported masses of quarks in 2011 give to the ratios m^ and rn^ an uncertainty of about 30% of their 
central value, therefore they were difficult to fit. The situation in 2012 has improved, in particular for the 
strange quark mass fhs the lattice determinations reduced its uncertainty to 5% of its central value. Since this 
is quite remarkable, we have also fitted the mass ratios using as limits the values obtained by considering just 
the theoretical determination of fhs, for which we obtain TOs(2 GeV) = 0.101 ± 0.011 GeV. 
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Variation of y^ with respect to m^ 
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Figure 2: The same as in Fig. [T] for x^ fit as a function of mc and irid- 
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Variation of y} with respect to m„ 
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Figure 3: The same as in Fig. [TJ except that now cos(/)i is allowed to vary in (—0.5, 0.5). 
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Variation of y^ with respect to m. 
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Figure 4: The same as in Fig. [31 except that now cos(/)i is allowed to vary in (—0.5, 0.5). 
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Variation of y} with respect to m„ 
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Figure 5: The same as in Figs.[T]|31 except that now coB(j)i is allowed to vary in the region (0.5, 1.0). 
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Variation of y^ with respect to m. 




0.0025 



0.003 



0.0035 



Variation of x ^^^^ respect to m. 




Figure 6: The same as in Figs.[T]|31 except that now cos(/)i is allowed to vary in the region (0.5, 1.0). 
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the theoretical determination of rris, we can see that the preferred region for the fit of fhc is 
quite similar to that of 2011. We can also see that the overall quality of the fit is better for the 
case of cos 01 if it is allowed to vary within (0.5, 1). In fact, if we allow cos0i to vary within 
(0, 1), the overall best fit is practically the same as that of the one when cos 01 is allowed to 
vary within (0.5, 1). 



Parameter 


Central value 


x' 


Values with restricted precision 


2 

X 


mu{Mz) 


1.72991 X 10"^ 




(1.73 ±0.75) X 10-5 




m, [Mz] 


3.46 X 10^3 




(3.46 ± 0.43) X 10-3 




md (Mz) 


1.12461 X 10~3 




(1.12 ±0.007) X 10-3 




rus [Mz] 


2.32 X 10-2 




(2.32 ± 0.84) X 10-2 




5u 


6.05040 X 10-2 




(6.05 ± 3.02) X 10-2 




Sd 


4.09162 X 10-2 




(4.09 ±2.59) X 10-2 




cos 01 


[Fixed] 


3.4 X 10-^ 




7.4 X 10-1 


mu{Mz) 


1.72960 X 10-^ 




(1.73 ±0.06) X 10-6 




m, (Mz) 


3.46008 X 10-3 




(3.46 ±0.31) X 10-3 




ma (Mz) 


9.19505 X 10-^ 




(9.20 ±0.72) X 10-4 




rris (Mz) 


2.08735 X 10-2 




(2.09 ±0.01) X 10-2 




K 


3.48158 X 10-2 




(3.48 ± 0.83) X 10-2 




Sd 


1.99291 X 10-2 




(1.99 ±0.63) X 10-2 




cos 01 


-1.42545 X 10-2 


1.32 X 10-5 


(-1.42 ±1.7) X 10-2 


2.4 X 10-2 


mu (Mz) 


1.71856 X 10-5 




(1.72 ±0.78) X 10-6 




rric (Mz) 


3.46176 X 10-3 




(3.46 ± 0.26) X 10-3 




rrid [Mz) 


1.05595 X 10-3 




(1.06 ±0.40) X 10-3 




rris {Mz) 


1.55660 X 10-2 




(1.56 ±0.72) X 10-2 




Su 


2.50428 X 10-2 




(2.50 ±6.18) X 10-2 




Sd 


4.09101 X 10-2 




(4.09 ± 7.04 X 10-2 




cos 01 


5.0 X 10-^ 


3.4 X 10-^ 


(5.0 ± 3.74) X 10-1 


1.6 X 10-1 



Table 8: Results of the fits for Case I, that is tlie case of a broken 5*3^ (E) Ssr symmetry. Note that when we 
restrict the precision of the fitted values, we observe a significant change in the value of x^. 



Interpretation of the parameters 5u and 6d- In this case, the symmetry breaking param- 
eter Zf, f = u,d, which measures the mixture of the singlet and the doublet representations of 
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Figure 7: \/5]; f ~ u blue (solid) curve, f ~ d orange (dashed) curve, as a function of Zj, for one of the 
three solutions of Eq. (P5|) . The plot on the left corresponds to the plot of 1999 from ref. jT^. while the plot 
on the right corresponds to taking the results, except the value of the parameters 5f, of the fit of case I for 
01 = 7r/2, Tab. [HI In this case, we have given as input the value of Zf and chosen the analogous solution to the 
1999 fit for Sf, using Eq. (|43|) . We see that both solutions are compatible, if we wanted to choose the value of 
the parameters Sf as a function of Zf. 



S3, was defined in ref. [T6l[T7] as 



'(Mf) 
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(Mf) 
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and it can be related to 6f tlirougli a cubic equation 



(42) 



F'Sf{Sf,Zf) 



^? 



Zf 



Zj + 1 



(nio 



mi. 



- f 
mi 



Zf{m{ 
m{Y 



m{){2 



0. 



V 



{l + 2{mi-m{))Zf)5) 



m^i. 



m,i 



mi){l 



m,i 



(43) 



Tlie equation {Zf + l)Fs^{6f, Zf) = is a linear function in Zf, hence given Sf there is only 
one solution for Zf. The x^ fits that we performed fitted the parameters 6f, hence fixing 
unequivocally the values of Zf. Had we chosen to fit Zf, we would have obtained three possible 
solutions for 6f. We find that ii Z = 0(10), the three solutions for 6f are roughly O(10~^), 
O(10~^) and 0(1), both ioi f = u and for f = d. In Fig. [7]we have plotted one of the solutions, 
that is close to 0(10""^), to Fs,{Sf,Zf) = 0, for two sets of data, that of 1999 [16] and the 
analogous solution for the case I of this study, when (pi is fixed to 7i/2. However, for this last 
set of data, we have used only the best fit points of the parameters rhi in Tab. [HI but not the 
value of 6f. Instead, once the parameters fhi are fixed, 5f{Zf) is computed as the solution of 
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Figure 8: On the left, Zf as a function oi 5f, f ~ u blue (solid) curve, f = d, orange (dashed) curve. The 
values obtained, from the fit for the case of 0i = 7r/2, for Su and Sd, 0.0605 and 0.0409 fix respectively particular 
values of Zf, 112.33 and 17.52. On the left, the parameter y/JJ as a function of Zf, f = u blue (solid) curve, 
f = d, orange (dashed) curve, for the data of 2012 for (pi = tt/2, corresponding to another solution of Eq. (P5|) . 
different from that of Fig. [H We have marked the points (17.52, VO.049) and (112.33, VO.0605). We can see 
that they correspond exactly to the solution of Eq. (|43[) , which is plotted here, since they lie along the line of 
F{6f, Zf) = 0, and not to the solution of Fig. [T] 



Eq. fj43|) that vanishes when Zf vanishes. In the first plot of Fig. [HI we show Zf as function of 
6f, computed from Eq. (H5]) for values of 6f close to the values obtained from the fit for case I 
with (pi fixed at 7r/2. The second plot in Fig. [H] shows ^y6f as a function of Zf. In these graphs 



6f{Zf) was chosen as the solution of Fs,{6f, Zf) = of O(10~^), in this way we check that it is 
indeed this solution the one determined in the x^ fit- 

In refs. [T6 | [T7 t l2"2 y 5 H l66] the solution to Eq. (H3|) was chosen such that Sf could represent a 
parameter of the breaking of 6*3, which would vanish in the limit of vanishing Zf. In the limit 
Z/ — )■ 0, Eq. (H3l) has two other solutions for 6f 



6f = l+ml and 6f = 1 - m(. 
Notice that neither of them satisfies the inequality 1 — m{ > 6f > |51) . 



(44) 



7.3.2 Cases II and III (three and four Higgs fields, respectively) 

In the cases of three and four Higgs fields, the functional form of the elements of the CKM 
matrix expressed as functions of the quark masses and the parameters /ig, 5u, and 6^ is the 
same. The only difference between the subcases A' and B' for the extended model with three 
Higgs fields (see Table 4) and among the subcases A', B', C, and D' for the extended model 
with four Higgs fields (see Table 5) is in the interpretation of the meaning of the parameters 
that occur as arguments in the elements of the CKM matrix in Eq. (137|) . which corresponds to 

and aj = mi + /Iq for /Iq 7^ 0. Thus, for this case, the relevant 
= 1, 2, 3, Su, Sd, and 02? which we denote as pi. Since the parameters 



taking Eq. ([MD with 0i 
parameters to fit are aj , 
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a/ are linear combinations of the mass ratios mj and the parameters JIq , we cannot fit JIq as 
part of the minimisation procedure, because then the fit would be underdetermined due to the 
linear dependence among its parameters. However, in order to give an interpretation to the 
parameters JIq, one possibility (i), is to let the parameters 5"/ to vary in the region 

rhi-3a,^f<a{<m{ + 3a^f, (45) 

and then check the compatibility of calculating /Iq such that 

/i^ = a"* - fhu, -^ Thc = a^- /Iq, (46) 

where •fhu, in the first expression of Eq. (l46l) . should lie within its experimental 3cr region. An 
analogous procedure is performed for the down-type quark sector. The values of a" and a'^ 
should correspond to the BFPs of the fit. Another possibility, (ii), is to fit the parameters pi, 
for a given value oi JIq 7^ 0, such that 

rhi - 3arh, + JJ^I < &u < rhi + Sa^^ + /i^, (47) 

with / = M for the up-type quarks and f = d for the down-type quarks. Note that the difference 
in the interpretations of /Xq for the cases above, Eqs. (l46|H7|) lies in the fact that for the case (i) 
there is no assumption on the value of /xq, but then /Iq necessarily has to be less than 6(7mu, 
while JIq has necessarily to be less than 60"^^- On the other hand, for (ii), in principle, there is 
no restriction on the value of the given JIq parameters. Although 5u and 5d have also a linear 
dependence on /Xq, Eq. (13T]) . since the definition involves another free parameter D^ ^ we can 
leave 5u and 5d to vary as completely free parameters. We performed two sets of fits, one using 
the 2012 values of the parameters rhi of Tab. [6] and the other the values of Tab. [TJ corresponding 
to the 2012 values of rhi when considering only the theoretical determination of rhg. We have 
allowed cos 02 to vary in the region (0, 1) and, we remind the reader that for this case the phase 
01, in Eq. (138|) . is equal to 0. The results of this fit are shown in Figs. l9]and[T0l 

In case (i), we can then calculate the value of /Iq from the experimental central value of m„, 
for which we have then JIq = —6 x 10~^ and consequently, from Eq. (I46p . fhc = 3.942 x 10~^, 
where we have used the values of au quoted in Tab. [9l Analogously for the d sector, we have 
fi'O — ^d ~ f^d = 3.8 X 10~^, as a consequence rhg = 2.04 x 10~^. Since we have the hierarchies 
rhu << fhc and m^ << rhg, while non-zero values of JIq and /Iq may be needed in this model to 
attain a best value of fhu and fhd-, respectively, concerning fhc and fhs, the impact is minimal. 

8 Conclusions 

We have studied the quark sector of different 6*3 models, either with one, three or four Higgs 
electroweak doublets. We presented the most general S'3-invariant Yukawa Lagrangian, which 
can describe these models. The structure of the Lagrangian gives rise to fermion mass matrices 
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Variation of y^ with respect to o„ 
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Figure 9: x^ ^s a function of (Tu and Cc, where coscfii is allowed to vary in the region (0.5, 1.0). The rest of 
the details are like those of Fig. [T] 
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Variation of y^ with respect to o^ 
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Figure 10: x^ a-s a function of ad and a-g, where cos0i is allowed to vary in the region (0.5, 1.0). The rest of 
the details are like those of Fig. [TJ 
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Parameter 


Central value 


x' 


Values with restricted precision 


x' 


Fit using the 2012 values of the parameters fhi 


^u{Mz) 
ac{Mz) 
ad{Mz) 
^s {Mz) 

Sd 

cos 02 


2.08977 X 10-6 
3.93180 X 10-3 
1.35949 X 10-3 
2.08443 X 10-2 
3.96726 X 10-2 
5.29260 X 10-2 
8.48776 X 10-1 


3.3 X 10-^ 


(2.09 ±0.19) X 10-6 
(3.93 ±0.007) X 10-3 
(1.36 ±0.004) X 10-3 
(2.08 ±0.02) X 10-2 
(3.97 ±0.35) X 10-2 
(5.29 ±041) X 10-2 
(8.49 ±0.22) X 10-1 


3.9 X 10-1 


Fit using the 2012 values o 


■ the parameters rhi (with mf-) 


^{Mz) 
^c{Mz) 
ad{Mz) 
as {Mz) 

Su 

5d 

cos 02 


2.17737 X 10-6 

3.94 X 10-3 
1.19392 X 10-3 
1.82432 X 10-2 
6.12747 X 10-2 
8.36979 X 10-2 
9.23028 X 10-1 


3.3 X 10-^ 


(2.18 ±0.35) X 10-6 
(3.94 ±0.007) X 10-3 
(1.19 ±0.009) X 10-3 
(1.82 ±0.02) X 10-3 
(6.13 ±0.41) X 10-2 
(8.37 ± 0.64 X 10-2 
(9.23 ±0.11) X 10-1 


7.3 X 10-2 



Table 9: Results of the fits for Cases II and III, that is the case of an SM invariant under an unbroken S'3 
symmetry. Note that when we restrict the precision of the fitted vahies, we observe a significant change in the 
value of x^ ■ 



of a generic form with a small number of free parameters, from which we were able to identify 
the conditions under which the two texture zeroes and Nearest Neighbour Interaction (NNI) 
mass matrices are obtained. In all cases we have provided exact, analytical formulas for the 
mixing angles of the CKM matrix in terms of quark mass ratios and a shift parameter //q. 
This line of work had been already developed in ref. [22|[73]. without referring to a particular 
model, where it was shown the usefulness of classifying mass matrix patterns according to their 
transformation properties under the group of permutations of three objects, S'3. There, it was 
also shown that a large class of phenomenologically succesful mass matrix forms are equivalent 
to two texture zeroes matrices. The reduction to these forms, two texture zeroes and NNI, and 
the fact that all CKM elements can be expressed as analytical relations in terms of quark mass 
ratios, allowed us to make a direct comparison of the models with the current experimental 
data. To this end, we performed a x^ fit of our theoretical expressions for the CKM mixing 
matrix to the experimentally determined values of the CKM matrix elements. 

In the case of the S'3 model with one Higgs electroweak doublet, which we have called S'3- 
SM, the S'3 symmetry has to be broken in order to give masses to all fermions. The resulting 
mass matrix, in a symmetry adapted basis, corresponds to a two zeroes texture. The value 
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of x^ of the fit to the CKM elements, is 1.6 x 10~^. In the case with three (Ss-SH) or four 
Higgs electroweak doublets the flavour symmetry is preserved. In these cases the resulting 
mass matrices correspond to either two zeroes textures or NNI ones, both known to be in good 
agreement with the phenomenology. The functional form of the CKM matrix elements is the 
same either with three or four Higgses, and the value of x^ is 3.9 x 10~^. 

It is worth noting that over the last decade there has been remarkable progress in reducing 
the uncertainties in the measurement of quark masses. Unlike one decade ago, presently it is 
no longer good enough to have a model that reproduces the hierarchy of fermion masses within 
the order of magnitude. At present, there are stringent limits on their values and so one must 
use statistical methods, such as a x^ fit, to measure the validity of a given model to reproduce 
the observed values for the CKM elements and the quark masses. The results of our x^ fits, 
show that the 6*3 models presented here reproduce with a remarkable accuracy the values of the 
CKM elements. The very good agreement between the 5*3 flavour symmetry models of quarks 
(presented in this work) and leptons [22l[5T] mixing, and the experimentally determined values 
of the corresponding mixing matrices, V^^^j^j and V^^f^^, gives a strong support to the idea 
that fermion masses and mixing might be related by a flavour permutational symmetry 5*3. 
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A Details of the rotation of mass matrices 

As we have mentioned in the main body of the text, the matrix Ai^^ makes explicit the 5*3 
transformations, however in order to diagonalise the mass matrix and extract the mixing matrix 
we perform a rotation and a shift as follows 

where A^^^jer' f*-*^ sub-cases A and A' of case III, it is explicitly given as 

/i{ + /i^c2(l -3t2) ^/sc(3 - t^) + ^^ 
-^L.^1 ^i{sc{?, - e) - ^ii ^{-^{c\l-?,e) ^i{/c |, (49) 
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while for the B and B' sub-cases of case III is given by 

/ /i( - /i{sc(3 - t^) _^/c2(l-3t2)+;xJ \ 

M^Hier=\ -fiic\l-3e)-fli fl{ + fiisc{3-t'') -fll/c , (50) 

V 4/c 4 J 

where c = cos0, s = sm6, and t = tan^. For sub-cases A, A', B, and B' of case II we just need 
to set /is = 0. Following Eq. ( TT6l) we can then identify that the condition for (1,3) and (3, 1) 
to vanish is 

tan6' = Wi/w2 (51) 

or 

tan^ = — W2/W1 (52) 

for sub-cases A and A' or B and B', respectively. 

The rotation in the Dirac fermion sector is unobservable, as long as we rotate both matrices, 
in the u and d or in the / and i^i sectors, by the same angle 6. Concerning the quark sector, 
the latter statement can be easily verified by diagonalising the matrices A^^ and the rotated 

matrices TZ{9) .^^ -^ Sz'^i^) fi2i 

ML,, = v-f^lLv-[^ = v-f7^(^)^,J7^(^).l2-MJ37^(^)^J7^(^).l2vi^ (53) 



It is then readily seen that the physical observables, contained in the CKM matrix, remain 
invariant 

v^cKM = v^n{e)l,,n{9),,X^. (54) 

Therefore, as long as we have the same rotation in both sectors, we preserve the matrix structure 
of the 5*3 symmetry. Now, we make a shift such that 

-^Ler = dhx3 + Mfmer, (55) 

where M.f Hier has the form of Eq. (123|) and hence we can proceed like in ref. [22] to diagonalise 
the mass matrix, which is explained in this work in Section 15.31 The matrix of Eq. (149|) can be 
identified with that of Eq. (12^ by assuming that 

r/ = z|r/|, Yi = Yi\ (56) 

where the first condition is needed such that the entries (1, 2) and (2, 1) correspond respectively 
to the complex conjugate of each other, and the second, such that the entries (2, 3) and (3, 2) 
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are also complex conjugate of each other. From Eqs. f l26|) and f H9|) we can see that the phase 
(pif is fixed by 

ta.n(j)if = —f = — 7 , (57) 

^ /i^sc(3-t2) |K/|w2Sc(3-t2)' 

or 

innSif = 7 = 7 , (58) 

^ -/i{c2(l-3t2) \Yi\wic'{l-?>t^) 

for sub-cases A and A' or B and B', respectively, which can be written in terms of the invariants 
of the matrix Ainier and the free parameter 5^ , which is the form we present in Eq. fl35l) . 

For sub-cases A' and B' of cases II and III, Tabs. ([T]) and (jl]), respectively, the form of the 
matrix can be reproduced just by assuming a rotation angle of ^ = 7r/6 oi 9 = 7r/3 for sub-case 
A' or B' , respectively, and without any necessity of imposing Hermiticity to the mass matrix. 

References 



[1 

[2 
[3 
[4 
[5 
[6 

[7; 

[8 

[9 

[10 

[11 
[12 

[13 



H. Fritzsch and Z.-z. Xing, Prog.Part.Nucl.Phys. 45, 1 (2000), [arXiv:hep-ph/9912358 

L. Velasco-Sevilla, J.Phys.Conf.Ser. 287, 012009 (2011). 

M. Gupta and G. Ahuja, Int.J.Mod.Phys. A26, 2973 (2011). 

H. Ishimori et al, Prog.Theor.Phys.Suppl. 183, 1 ^2010). larXiv: 1003.35521 

G. Altarelh and F. Feruglio, Rev.Mod.Phys. 82, 2701 f2010). laiOGv: 1002.02111 

M. Hirsch et al, (2012). E?Xrv:1201.5525[ 

R. Gatto, G. Sartori, and M. Tonin, Phys.Lett. B28, 128 (1968). 

N. Cabibbo and L. Maiani, Phys.Lett. B28, 131 (1968). 

H. Pagels, Phys.Rev. Dll, 1213 (1975). 

S. Weinberg, Trans.New York Acad.Sci. 38, 185 (1977). 

F. Wilczek and A. Zee, Phys.Lett. B70, 418 (1977). 

H. Fritzsch, Phys.Lett. B70, 436 (1977). 

A. Ebrahim, Phys.Lett. B73, 181 (1978). 

39 



[14] R. N. Mohapatra and G. Senjanovic, Phys.Lett. B73, 176 (1978). 

[15] H. Fritzsch, Phys.Lett. B73, 317 (1978). 

[16] A. Mondragon and E. Rodriguez- Jauregui, Phys.Rev. D59, 093009 (1999), 



arXiv:hep-ph/9807214 



[17] A. Mondragon and E. Rodriguez- Jauregui, Phys.Rev. D61, 113002 (2000), 



arXiv:hep-ph/9906429 



[18] S. Morisi, (2006), arXiv:hep-ph/0604106 



[19] F. Ferugho and Y. Lin, NucLPhys. B800, 77 (2008), arXiv:0712. 15281 

[20] T. Kobayashi, Y. Omura, and K. Yoshioka, Phys.Rev. D78, 115006 (2008), 
larXiv: 0809. 30641 

[21] R. Jora, J. Sc hechter, and M. Naeeni Shahid, Phys.Rev. D80, 093007 (2009), 
larXiv:0909.4414[ 

[22] J. Barranco, F. Gonzalez Canales, and A. Mondragon, Phys.Rev. D82, 073010 (2010), 
larXiv: 1004.37811 

[23] Z.-z. Xing, D. Yang, and S. Zhou, Phys.Lett. B690, 304 (2010), arXiv: 1004.42341 

[24] S. Zhou, Phys.Lett. B704, 291 r20in. la7xrv:1106.4808[ 

[25] D. Meloni, JHEP 1205, 124 ^2012). l^?Xrv:1203.3126[ 

[26] S. Dev, R. R. Gautam, and L. Singh, Phys.Lett. B708, 284 (20 12) . li^XTv: 1201. 37551 

[27] H. Benaoum, (;2013). laiOGv: 1302.09501 

[28] S. Pakvasa and H. Sugawara, Phys.Lett. B73, 61 (1978). 

[29] E. Dernian, Phys.Rev. D19, 317 (1979). 

[30] D. Wyler, Phys.Rev. D19, 330 (1979). 

[31] J.-M. Frere, Phys.Lett. B80, 369 (1979). 

[32] R. Yahalom, Phys.Rev. D29, 536 (1984). 

[33] E. Ma, Phys.Rev. D43, 2761 (1991). 



[34] L. J. Hall and H. Murayama, Phys.Rev.Lett. 75, 3985 (1995), [a^OCrv:hep-ph/9508296 

40 



[35] L. Lavoura, Phys.Rev. D61, 077303 (2000), [a^Xw:hep-ph/9907538 



[36] Y. Koide, Phys.Rev. D60, 077301 (1999), ]arXiv:hep-ph/9905416 Revised version. 



[37] J. Kubo, A. Mondragon, M. Mondragon, and E. Rodriguez- Jauregui, Prog.Theor.Phys. 
109, 795 (2003), [a^Xw:hep-ph/0302196J 



[38] J. Kubo, H. Okada, and F. Sakamaki, Phys.Rev. D70, 036007 (2004), 



arXiv:hep-ph/0402089 



[39] S.-L. Chen, M. Frigerio, and E. Ma, Phys.Rev. D70, 073008 (2004) , [^?Xw:hep-ph/0404084l 



[40] Y. Koide, Phys.Rev. D73, 057901 (2006), arXiv:hep-ph/0509214 



[41] T. Kiniura, Prog.Theor.Phys. 114, 329 (2005). 

[42] T. Araki, J. Kubo, and E. A. Paschos, Eur.Phys.J. C45, 465 (2006), 



arXiv:hep-ph/0502164 



[43] A. Mondragon, M. Mondragon, and E. Peinado, J.Phys. A41, 304035 (2008), 
larXiv:0712. 17991 

[44] S. Kaneko, H. Sawanaka, T. Shingai, M. Tanimoto, and K. Yoshioka, (2007), 
arXiv:hep-ph/0703250 , 

[45] A. Mondragon, M. Mondragon, and E. Peinado, AIP Conf.Proc. 1026, 164 (2008), 
[i7xrv:0712.2488. 

[46] A. Mondragon, M. Mondragon, and E. Peinado, Phys.Rev. D76, 076003 (2007), 
[a?Xrv:0706.0354i 

[47] O. F. Beltran, M. Mondragon, and E. Rodriguez- Jauregui, J.Phys.Conf.Ser. 171, 012028 
(2009). 

[48] G. Bhattacharyya, P. Leser, and H. Pas, Phys.Rev. D83, 011701 (2011),[arXml0Q63597| 

[49] T. Teshima and Y. Okumura, Phys.Rev. D84, 016003 (;2011). [a?Xrv:1103.6127[ 

[50] T. Teshima, Phys.Rev. D85, 105013 (20 12) . [^ilGv: 1202.45281 

[51] F. Gonzalez Canales, A. Mondragon, and M. Mondragon, Fortschritte der Physik (2012), 
larXiv: 1205.47551 

[52] Particle Data Group, J. Beringer et al, Phys. Rev. D86, 010001 (2012). 

[53] G. Fogh et al, Phys.Rev. D86, 013012 (2012). liDCTv: 1205.52541 

41 



[54 
[55 
[56 
[57; 
[58 
[59 
[60 

[61 
[62 
[63 
[64 
[65 
[66 

[67; 
[68 

[69 

[7o; 

[71 
[72 



D. Forero, M. Tortola, and J. Valle, Phys.Rev. D86, 073012 (2012), arXiv: 1205 .40781 

Particle Data Group, K. Nakamura et a/., J.Phys.G G37, 075021 (2010). 

A. Mondragon et al, Work in progress. 

J. Girrbach, p. 83 (;2012V [i7xrv:1208.5630[ 

S. F. King and C. Luhn, (2013). laTXTvilSOl. 13401 

D. Emmanuel-Costa, C. Simoes, and M. Tortola, (;2013). farXiv:1303.56991 

J. Gomez-Izquierdo, F. Gonzalez Canales, and M. Mondragon, A Grand Unified Model 
with Q6 as the Flavour Symmetry. To appear in the proceedings of PASCOS 2012. 

CMS Collaboration, S. Chatrchyan et al, (;2013). HarXiv:1302.1764[ 



N. Haba and K. Yoshioka, Nucl.Phys. B739, 254 (2006), |arXiv:hep-ph/0511108 



ATLAS Collaboration, G. Aad et al, Phys.Lett. B716, 1 (;2012V E?Xrv:1207. 72141 
CMS Collaboration, S. Chatrchyan et al, Phys.Lett. B716, 30 f2012). EFXTv: 1207. 72351 
CMS Collaboration, S. Chatrchyan et al, (;2013). l [a?Xrv:1303.4571[ 
A. Mondragon and E. Rodriguez- Jauregui, Rev.Mex.Fis. 46, 5 (2000), 



arXiv:hep-ph/0003104 , 



K. Harayama and N. Okamura, Phys.Lett. B387, 614 (1996), arXiv:hep-ph/9605215 



K. Harayama, N. Okamura, A. Sanda, and Z.-Z. Xing, Prog.Theor.Phys. 97, 781 (1997), 



arXiv:hep-ph/9607461| 



F. Gonzalez-Canales, A. Mondragon, U. Saldafia Salazar, and L. Velasco-Sevilla, (2012), 
larXiv: 1210.02881 

K. Chetyrkin, J. H. Kuhn, and M. Steinhauser, Comput.Phys.Commun. 133, 43 (2000), 



arXiv:hep-ph/0004189 



Z.-z. Xing, H. Zhang, and S. Zhou, Phys.Rev. D77, 113016 f 2008 V [a7xrv:0712. 14191 

R. Brun and F. Rademakers, Proceedings AIHENP'96 Workshop, Lausanne, Sep. 1996, 
Nucl. Inst. & Meth. in Phys. Res. A , 81 (389 (1997)). 



[73] F. Gonzalez Canales and A. Mondragon, J.Phys.Conf.Ser. 287, 012015 (2011), 
larXiv:1101.3807[ Presented at XIV Mexican School on Particles and Fields, 4-13 November 
2010, Moreha Mexico. 



42 



